function pds=pd_plot(equity_mat,optimal_coupon,n_p0,setts)
% Plot default thresholds for p=1,...,n_p0

    pds = zeros(n_p0+1,1);
    for k=1:n_p0+1
        pds(k)=sum(equity_mat(find(round(setts.C,4)==round(optimal_coupon(k),4)),:)==0);
    end

end